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In this paper on the basis of the generalized f(R) gravity model with arbitrary coupling between 
geometry and matter, four classes of f(R) gravity models with non-minimal coupling between geom- 
etry and matter will be studied. By means of conditions of power-law expansion and the equation 
of state of matter less than —1/3, the relationship among p, w and n, the conditions and the candi- 
date for late-time cosmic accelerated expansion will be discussed in the four classes of f(R) gravity 
models with non-minimal coupling. Furthermore, in order to keep considering models to be realistic 
D ■ ones, the DolgovCKawasaki instability will be investigated in each of them. 
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Introduction 



According to recent observational data sets [1-3], our current universe is flat and undergoing a phase of the 
accelerated expansion which started about five billion years ago. To explain this phenomena, a variety of models have 
00 ■ been proposed which may be divided into three broad classes. First, it is possible that there is some undiscovered 
• property in our existing model of gravity and matter that leads to acceleration at the present time. In thses scenarios, 



one might consider including the existence of a tiny cosmological constant and the possibility of the backreaction of 
cosmological perturbations. 

Second is the idea that the Universe is dominated by an exotic component with large negative pressure, usually 
, referred to as dark energy. The simplest form of dark energy is cosmological constant A which would encounter 
fine-tuning problem and coincidence problem. Other valid dark energy models are provided by scalar fields, such as: 
■ Quintessence [4. 5], which is introduced to solve the coincidence problem and characterized by the equation of state 
(EOS) Wde between -1 and -1/3 (namely, —1 < Wd e < —1/3); Phantom (ghost) field[6], which owns a negative kinetic 
energy and characterized by the EOS Wde less than -1 (namely, Wde < — 1); Tachyon field[7, 8]which can act as a 
source of dark energy depending upon the form of the tachyon potential, and so on. Other scenarios on dark energy 
include brane world[9], generalized Chaplygin gas[10], holographic dark cnergy[ll], etc. Unfortunately, up to now a 
satisfactory answer to the question that what dark energy is and where it came from has not yet to be obtained. 
Finally, eliminating the need of dark energy, one may consider modified theories of gravity in which the late- 
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time cosmic accelerated expansion can be realized by an infrared modification. There are numerous ways to deviate 
from Einstein's theory of gravity. The most well-known alternative to General Relativity (GR) is scalar-tensor 
theory [12, 13]. There are still numerous proposals for modified theories of gravity in contemporary literature, such as 
DGP (Dvali- Gabadadze-Porrati) gravity [14], braneworld gravity [15], TeVeS (Tensor- Vector- Scalar) [16], f(R) theories 
of gravity [17], Einstein- Aether theory [18] and so on. 

Among these theories, f(R) gravity is very competitive. Here f(R) is an arbitrary function of the Ricci scalar R. 
One can add any form of R in it, such as l/i?[19] (the simplest one), lni?[20], positive and negative powers of i?[21], 
Gauss-Bonnet invariant [22], etc. The more general forms of f(R) can be considered including coupling between f(R) 
and scalar [23], multidimensional 1/R theory [24] and so on. In f(R) theories of gravity, the expansion history of the 
universe is naturally explained by the fact that some gravitational terms which support the inflation at early-time 
universe, while other terms which cause the cosmic acceleration at late-time universe. It is worth stressing that 
considering some additional conditions, the early-time inflation and late-time acceleration can be unified by different 
role of gravitational terms relevant at small and at large curvature. However, f(R) gravity is not perfect because 
of containing a number of instabilities. For instance, the theory with 1/R may develop the instability [25]. But by 
adding a term of R 2 to this specific f(R) model, one can remove this instability[20, 21]. For more general forms of 
f(R), the stability condition /" > can be used to test f(R) gravity models [26]. 

Recently, a general model of f(R) gravity has been proposed in Ref.[27], which contains a non-minimal coupling 
between geometry and matter. This coupling term can be considered as a gravitational source to explain the current 
acceleration of the universe. The viability criteria for such a theory was recently discussed in Refs.[27, 28]. However, 
a more general model, in which the coupling style is arbitrary and the Lagrangian density of matter only appears 
in coupling term, has been proposed in Rcf.[29] and it can represent the former case. The purpose of this paper 
is to discuss the conditions for latc-time behaviour in f(R) gravity with non- minimal coupling between matter and 
geometry. 

This paper is organized as follows. In the next section, for the general f(R) gravity models with arbitrary and 
non-minimal couplings, the field equations, the energy conditions and the DolgovCKawasaki instability will be given, 
respectively. In this class of models, the energy-momentum tensor of matter is generally not conserved due to the 
appearance of an extra force as mentioned in Rcf. [27]. The conditions for late-time cosmic accelerated expansion and 
the instability in f(R) gravity with non-minimal coupling will be discussed in sections 3 and 4, respectively. Four 
classes of models will be taken into consideration in those two sections. Summary is given in the last section. 

2. THE GENERAL f(R) GRAVITY WITH COUPLING BETWEEN MATTER AND GEOMETRY 

A more general action in f(R) gravity, in which the coupling style between matter and geometry is arbitrary 
and the Lagrangian density of matter only appears in coupling term, is given by 



where we have chosen k = 8irG = c = 1, which we shall adopt hereafter. fi{R) (i = 1,2) and G(L m ) are arbitrary 
functions of the Ricci scalar R and the Lagrangian density of matter respectively. When f2(R) = 1 and G(L m ) = L m , 
we obtain the general form of f(R) gravity with non-coupling between matter and geometry. Furthermore, by setting 




(1) 
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fi(R) = R, action (1) can be reduced to the standard General Relativity (GR). 
Varying the action (1) with respect to the metric g^ yields the field equations 

FiORj-R^ - yi(R)g„ v + {g^a - V^^F^R) = -2G(L m )F 2 (R)R^ 

-2(g MV D - V^ v )G{L m )F 2 {R) - f 2 (R)[K(L m )L m - G(L m )]g^ (2) 

+f 2 (R)K(L m )T fiV , 

where □ = g^V^Vv, Fj(R) = dfi(R)/dR (i = 1,2) and K(L m ) = dG(L m )/dL m . The energy-momentum tensor of 
matter is defined as: 



2 S(^L m ) 
~~g Sg^ 



Assuming that the Lagrangian density of matter L m only depends on the metric tensor components and not on its 
derivatives, we obtain 

T^ v = L m g^ v - 2^^. (4) 

The trace of the field equations (2) reads 

3D[*i(JZ) + 2G(L m )F 2 (R)} + [FAR) + 2G(L m )F 2 (R)]R 

( 5 ) 

-2h(R) + if 2 (R)[K(L m )L m - G(L m )\ = K(L m )f 2 (R)T, 

where T = Ttf. 

By taking the covariant divergence of Eq.(2) and using the mathematical identity V M [f'(R)R flv — hf(R)g flv + 
(g^U - V^V^/GR)] = [29], here f'(R) = df/dR, we have 

VT^ = 2^ln[f 2 (R)K(L m )]-^, (6) 

from which we see that the conservation of the energy-momentum tensor of matter is violated due to the coupling 
between matter and geometry. However, once the L m is given, by choosing appropriate forms of G(L m ) and f 2 (R), 
one can construct, at least in principle, conservative model with arbitrary matter-geometry coupling. 

In order to keep the energy density is positive and cannot flow faster than light, the generalized energy conditions, 
namely, the strong energy condition (SEC), the null energy condition (NEC), the weak energy condition (WEC) and 
the dominant energy condition (DEC), should be taken into consideration, which forms can be derived as follows (see 
Rcf. [30] for more details): 

P + 3p-jfclfi- (fl + *Gft)R] + Z^{HR + R) 
+3j^R 2 + 6j^gt(G'' L m f 2 + L m G' f 2 + 2/2' RG'L m 
+ti"R 2 G + f 2 RG) + %-M^(G'L m f 2 + f 2 RG) 
+£(G'L m -G) >0, (SEC) 



p + p+(HR + R)^ + j^R? + jfa{G»L m 2 ft+ 

L m G'f 2 + 2f±RG'L m + f 2 "R 2 G + f 2 RG) - (8) 
^(G'L m f 2 + ti'RG) > 0, (NEC) 
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P-P+jZplfl ~ (/( + 2G/^)i?] - (5HR + R) 7 g 7 - 
JTg 7 ^ 2 ~ j^g 7 ^" L m / 2 + L m G f 2 + 2/ 2 RG'L m + 
f 2 'R 2 G + f 2 RG) - jP§r(G r L m f2 + f 2 RG) - J7 
(G'L m - G) > 0, {DEC) 



{G'Lrnf'2 + flRG) - £r (G'L m —G)>0, (WEC) 

where the dot denotes differentiation with respect to cosmic time. Moreover, by using condition f"(R) + 
2G(L m )f 2 (R) > 0, one can test the Dolgov-Kawasaki instabilities for this class of models. 

When G(L m ) = L m and rescales the function f2(R) as 1 + Xf2(R), the action (1) and the field equations (2) can 
be changed into 

S=f {\h(R) + [1 + \f2{R)]L m }V~gd i x 1 (11) 

- \h{R)g^ + - V M V„).Fi(iZ) = -2XF 2 (R)L m R^ 

+2A(V M V„ - 5 ^D)L m F 2 (i?) + [1 + Xf 2 (R)]T^, 

Above expressions are just the action and the field equations in f(R) gravity with non-minimal coupling between 
geometry and matter. Moreover, by means of the generalized Bianchi identities V M G M „ = (here, G M „ is the Einstein 
tensor), Eq.(6) can be given as: 

XFn 

= % [g^L m - T^Jv^iZ. (13) 

I + A/2 

It follows that the non-minimal coupling term results in a non-trivial exchange of energy and momentum between 
geometry and matter[31, 32]. Note that according to Eq.(13), the conservation of the energy-momentum tensor can 
be verified if fziR) is a constant or the Lagrangian density of matter is not an explicit function of the metric. 



3. THE CONDITIONS FOR LATE-TIME COSMIC ACCELERATED EXPANSION IN f(R) GRAVITY 

WITH NON-MINIMAL COUPLING 

In the following, we focus on the conditions for late-time cosmic accelerated expansion in f(R) gravity with 
non-minimal coupling between geometry and matter. In this section, the form of the action is taken to be Eq.(ll) 
and, for simplicity, we consider L m is opposite to the energy density of perfect fluid[32], i.e., 

L m = -p=-p a-^ 1+w \ (14) 

where w is the equation of state of perfect fluid and is assumed to be a constant. The energy- momentum tensor is 
taken as: 



Tpv = {p + pjUfjUv + pg^ v , 



(15) 
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where p and p denote the energy density and the pressure respectively. The form of the FRW metric is chosen as 

ds 2 = -dt 2 + a 2 (t)dX 2 , (16) 

where a(t) is the scale factor and dX 2 contains the spacial part of the metric. Using this metric, we can obtain 
R = 6(2H 2 + H), where H = a(t)/a(t) is the Hubble expansion parameter. 

It is known that under the conditions cither power-law expansion or the equation of state of matter less than 
— 1/3, late-time cosmic accelerated expansion occurs. To exemplify how to use these conditions to realize the phase 
of accelerating expansion in f(R) gravity with non- minimal coupling, Firstly, we concentrate on two simple classes of 
models. 

(1) Let 

fi(R) = R, / 2 (i?) = -AR- n + BR 2 , (17) 

where A and B are arbitrary constants. Then, Eq.(12) becomes into 

3H 2 = -p a" 3(1+, " ) [l + 6\{H 2 + H){2ABH 2 + 12BH+ 
An(12H 2 + 6H)- 1 - n ) + \(36B(2H 2 + H) 2 - A(12H 2 + 6H)~ n )]. 



We assume the solution of Eq.(18) is a = aot p and then have H = 2, H = — J?. Substituting these relations into 

2 

3(l+w) 



Eq.(18), we find there are three kinds of possible relationships among p, w and n, namely, p = 3("+ui) ' P ~ — ~ 



and p = 3 /^|_ M) \ ■ Under the condition of power- law expansion (i.e., p > 1), the corresponding regions of w are 
w < 2 ^" 3 f 1 - ) — 1 for p — §fep^ ; w < — 1/3 for p = 3 ^ w ^ and w < —5/3 for p = 3 /^ w \ , respectively. Furthermore, by 
considering the equation of state of matter less than —1/3 (i.e., w < —1/3), we can obtain that when p = , the 



range of parameter n is n < and n/-l. It is easy to see that there is no constraint on n when cases p = „ (1 2 ■. and 



3(l+tu) 1 
2 

3(l+w) 



P = 3{i+w) ■ ^ 0T ^ ne case P = 3,{i+ w) ' * ne en?ec tive quintessence regime (—1 < w < —1/3) emerges when — 1 < n < 
and the effective phantom regime (w < —1) emerges when n < — 1. The candidate for late-time cosmic accelerated 
expansion can be either the effective quintessence or the effective phantom, when p = 3^^^ • 

(2) Another choice for functions fi(R) and f2(R) are 

fi(R) = R, f2(R) = c J n + 1 , (19) 
where c\ and C2 are constants. Then the FRW equation is changed into 

m 2 = - P0 a-^[i + 6 " ClwA(g2 + g)(2g2 t g) ' 1+ " + g^gg^+jr 

[1 + 6 n c 2 (2H 2 + H) n } 2 1 + 6' 1 c 2 (2H 2 + H) n 

By calculations and analysis, the relationship among p, w and n, condition and candidate for late-time cosmic 
accelerated expansion are shown in Table I. 

It is worth stressing that the de Sitter stage is impossible in both models because there is a scale factor a in the 
FRW equation. If the function f2(R) vanishes, the standard FRW equation would be reproduced. Above forms of 
fi(R) have been discussed in Refs.[21, 33]. Next, we focus on other two complicated models. 

(3) Following Rcf. [21], let us take the following explicit choice for functions fi(R) and f2{R) as: 

f 1 (R)=R-AR- n + BR 2 , h{R) = —AR~ n + BR 2 , (21) 
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Relationship 


Condition 


Candidate 


The effective quintessence 


The effective phantom 


_ _ 2(i-«) 

P 3(l+to) 


n > and n / 1 


< n < 1 


n > 1 


2(l-2n) 
P 3(l+io) 


n > and n 1/2 


< n < -1/2 


n > 1/2 


P — 3(l+«l) 


id < -1/3 


All n 


All n 



TABLE I: The relationship among p, w and n, condition and candidate for late-time cosmic accelerated expansion in case 
fi{R) = R, h{R) = TTrSxt- 



where A and B are arbitrary constants. Then the equation (12) can be expressed as: 

\{\2H 2 + 6H + 36B(2H 2 + Hf - A(12H 2 + 6H)- n - 6(H 2 +H)[1 + 24BH 2 + 

12BH + An{l2H 2 + 6H)- 1 -"}} = -p a-^ 1+w ^{l + 6X(H 2 + H)[2ABH 2 + 12BH+ (22) 

An(l2H 2 + 6H)- 1 -"} + X[36B(2H 2 + H) 2 - A(12H 2 + 6H)~ n ]}. 

By the same method as above, we find there are four kinds of possible relationships among p, w and n, i.e., 
p = 3 ( 1 ^_ tu ) ; P = — 4^3^ i P = 3(i4-ti,) an d P = — 2n. Under conditions of power-law expansion and the equation of state 
of matter less than —1/3, above results are turned into w < 1/3 for p = 3 ^ w ^ , n > —3/2 for p = — 4 ^ w , w < —1/3 
for p = arL d n < — 1/2 for p = — 2n, respectively. It is clear that the candidate for late-time cosmic accelerated 

expansion is among dust, the effective quintessence and the effective phantom when p = 3 ( 1+w ) ■ But except dust 
when p = -rrr. — r. Note that there is no constraint on n in both cases. For the case p = — , 2r t , either the effective 
quintessence regime emerges when —3/2 < n < —1/2 or the effective phantom regime emerges when n > —1/2 and 
n =/= 0. Obviously, the late-time cosmic accelerated expansion is independent of matter when p = —2n. 
(4) Another choice for functions fi(R) and f2(R) are 

f 1 (R) = R+—±- — , h{R) = — 1 - , (23) 

where c\ and c 2 are constants. In this case, the FRW equation can be given as: 



■W2H 2 + 6H+ 6" Cl (2g 2 +g)" _ 6(H 2 , , 6- 1 +" Cl »(2g 2 +g)- 



„ „-3(l+«;)ri , 6 n c 1 n\(H 2 +H)(2H"+H)- 1 + n 6" Cl A(2g 2 + g)" i 
W !■ [l+6"c 2 (2ff 2 +H)«] 2 l+6"c 2 (2ff 2 +_f/)" J ' 



(24) 



The corresponding relationship among p, w; and n, condition and candidate for late-time cosmic accelerated expan- 
sion are shown in Table II. 

Obviously, the de Sitter stage is also impossible in both models and the reason is as the same as the ones in above 
two simple cases. If the function f2(R) vanishes, the modified gravity with non-coupling can be reproduced. Above 
forms of fi(R) have been discussed in Rcfs.[21, 34]. 

From the above discussions, it is easy to see that the results in complicated models are more interesting than the 
simple ones. 

For the four models mentioned above, the transition from matter- dominated phase to the acceleration phase 
discussions without non-minimal coupling have been made in Ref.[35]could be realized as follows. Since the Hubble 
parameter can be expressed as H = p/t, Ricci scalar R turns into R = 6p(2p — l)/t 2 . If < p < 1, the early universe 
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Relationship 


Condition 


Candidate 


The effective quintessence 


The effective phantom 


a{i-\-W) 


w < -1/3 


All n 


All n 


4n + 2 
F 3{1+uj) 


n < and n / -1/2 


-1/2 < n < 


n < -1/2 


„ _ 2(n + l) 
r ~ 3(1+™) 


n < and n/-l 


-1 < n < 


n < -1 


f 3(1+™) 


n < 1 and n/0 


< n < 1 


n< 


n — 4n 
r ~ 3(1+™) 


n < 1/2 and n / 


< n < 1/2 


n < 



TABLE II: The relationship among p, and n, condition and candidate for late-time cosmic accelerated expansion in case 



is in deceleration phase, which corresponds to matter-dominated phase with p = 2/3, and if p > 1, the late universe 
is in acceleration phase. 

4. THE INSTABILITY OF f(R) GRAVITY WITH NON-MINIMAL COUPLING 

A viable modified gravity model must pass Newton law, solar system test and instability conditions [21, 36]. 
There are in principle several kinds of instabilities to considcr[37]. DolgovCKawasaki instability [38] is one of them. 
Below, we will focus on this instability. According to Ref.[30], the DolgovCKawasaki criterion in f(R) gravity with 
non-minimal coupling between matter and geometry is 

f['(R)+2\L m f^(R)>Q (25) 

where A is a constant and L m is the Lagrangian density of matter. For simplicity assuming A, B, C\, C*2 are positive 
constants, the DolgovCKawasaki criterions for above four discussed models are as follows: 

n < or n < — 1, for model 1 and model 3, (26) 



1 -|- C R n 

n < — , for model 2 and model 4 (27) 

~ 1 - C 2 R n 

By means of analysis of model 1, model 3 and Eq.(26), we find that they could be realistic candidates for late- 
time cosmic accelerated expansion without DolgovCKawasaki instability. By taking Table 1, Table II and Eq.(27) into 
consideration, model 2 and model 4 would be realistic candidates if 1 + C2-R n /1 — CiR n > and 1+C2-R"/1— C2.R" < 1 
respectively Otherwise there is no interesting in model 2. Furthermore, the DolgovCKawasaki instability will be 
emerged in model 4, if 1 + C 2 R n /l - C 2 R n < 1/2 (i.e. 1/2 < n < 1) for case p = 2n/(l + w), < n < 1 for case 
p = 2n/ (1 + w) and < n < 1/2 for case p = An/ (1 + w), respectively. 



5. SUMMARY 



Up to now, we have discussed the conditions for late-time cosmic accelerated expansion and the DolgovCKawasaki 
instability in f(R) gravity with non-minimal coupling between geometry and matter. For simplicity, we chose the 
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form of the Lagrangian density of matter as opposite to the energy density of perfect fluid. The relationship among 
p, w and n has been given in each class of models. By using the conditions of power-law accelerated expansion, 
the equation of state of matter less than —1/3 and the DolgovCKawasaki criterion, the range of the parameter n is 
concretely constrained. Either the effective quintessence regime or the effective phantom regime would emerge by 
choosing n properly. It is easy to see that the results in complicated models are more interesting than the simple ones. 
It is demonstrated that the de Sitter stage would not realize in all considering models because there is a scale factor 
a in FRW equation. Essentially, this is due to the special choice of the Lagrangian density of matter. Other forms of 
the Lagrangian density of matter could be considered in the similar fashion of non-minimal gravitational coupling. 



1] Adam G. Riess et al, Astron. J. 116, 10091998 
2] S. Perlmutter et al., Astrophys. J. 517, 565 1999 
3] M. Tegmark et al, Astrophys. Rev. D 69, 103501 (2004) 

M. Tegmark et al, Astrophys. J. 606, 702 (2004) 
4] Wetterich C, Nucl. Phys. B 302, 668 (1988) 
5] Amendola L., Pfrys.Rev. D 62, 043511 (2000) 
3] Caldwell R. R., Phys. Lett. B 545, 23-29 (2002) 
7] Padmanabhan T., Phys. Rev. D 66, 021301 (2002) 

S] Bagla J. S., Jassal H. K., Padmanabhan T., Phys. Rev. D 67, 063504 (2003) 
9] Randall L., Sundrum R., Phys. Rev. Lett. 83, 4690 (1999) 
3] Bento M. O, Bertolami O., Sen A. A., Phys. Rev. D 66, 043507 (2002) 
1] Li M., Phys. Lett. B 603, 1 (2004) 
2] C. Brans, R.H. Dicke, Phys. Rev. 124, 925 (1961) 

3] V. Faraoni, Cosmology in Scalar- Tensor Gravity (Kluwer Academic, Dordrecht, 2004) 
4] Dvali et al., Phys. Lett. B 485, 208-214 (2000) 
5] R. Maartens, Living Rev. Rel 7, 7 (2004) 
S] J.D. Bekenstein, Phys. Rev.D 70, 083509 (2004) 

7] S. Nojiri, S.D. Odintsov, Int.J.Geom.Meth.Mod.Phys. 4, 115-146 (2007) 
S] T. Jacobson, D. Mattingly, Phys. Rev.D 64, 024028 (2001) 
9] S. Capozziello, Int. J. Mod. Phys.D 11, 483 (2002) 
3] S. Nojiri, S.D. Odintsov, GRG 36, 1765 (2004) 
1] S. Nojiri, S.D. Odintsov, Phys. Rev.D 68, 123512(2003) 

Gianluca Allemandi, Andrzej Borowiec, Mauro Francaviglia, Sergei D. Odintsov, Phys. Rev.D 72, 063505 (2005) 
T. Inagaki, S. Nojiri, S. D. Odintsov, JCAP 0506, 010 (2005) 
M.C.B. Abdalla, S. Nojiri, S. D. Odintsov, Class.Quant.Grav. 22, L35 (2008) 
[22] S. Nojiri, S.D. Odintsov, Phys.Lett.B 631, 1-6 (2005) 
[23] S. Nojiri and S.D. Odintsov, arXiv:astro-ph/0403622 

G. Cognola, S. Zerbini, arXiv:gr-qc/0407103 
[24] U. Guenther, A. Zhuk, V Bezerra, C. Romero, Class.Quant.Grav. 22, 3135-3167 (2005) 
[25] A.D. Dolgov, M. Kawasaki, Phys.Lett. B 573, 1-4 (2003) 
[26] V. Faraoni, Phys. Rev.D 74, 104017 (2006) 



9 

[27] O. Bertolami, C.G. Boehmer, T. Harko, F.S.N. Lobo, Phys. Rev.D 75, 104016 (2007) 

[28] V. Faraoni, arXiv:0710.1291 [gr-qc] 

[29] T. Harko, Phys.Lett.B 669, 376-379 (2008) 

[30] J. Wang, Y. B. Wu, Y. X. Guo, W. Q. Yang, L. Wang, Phys. Lett. B 689, 133-138 (2010) 
[31] O. Bertolami, J. Pramos, Class. Quant. Grav. 25, 245017 (2008) 

[32] O. Bertolami, T. Harko, Francisco S. N. Lobo, Jorge Pramos, arXiv:0811.2876 [gr-qc] 

J.D. Brown, Class. Quant. Grav. 10, 1579 (1993) 
[33] W. Hu, I. Sawicki, Phys.Rev.D 76, 064004 (2007) 

[34] S. Nojiri, S. D. Odintsov, P. V. Tretyakov, Prog.Theor.Phys.Suppl. 172, 81-89 (2008) 
[35] S. Nojiri, S. D. Odintsov, Phys.Rev.D 74, 086005 (2006) 

K. Bamba, S. Nojiri, S. D. Odintsov, JCAP 0810, 045 (2008) 
[36] S. Nojiri, S. D. Odintsov, Phys.Lett.B 652, 343-348 (2007) 

S. Nojiri, S. D. Odintsov, Phys.Lett.B 657, 238-245 (2007) 

S. Nojiri, S. D. Odintsov, Phys.Rev.D 77, 026007 (2008) 
[37] Alvaro Nunez, Slava Solganik, arXiv:hep-th/0403159 
[38] A. D. Dolgov, M. Kawasaki, Phys. Lett. B 573 (2003) 1 



